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Density profiles and substructure of dark matter halos: 
converging results at ultra-high numerical resolution. 
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(N . ABSTRACT 

f^ i- Can dissipationless N-body simulations be used to reliably determine the structural and sub- 

l' ' structure properties of dark matter halos? A large simulation of a galaxy cluster in a cold 

dark matter universe is used to increase the force and mass resolution of current "high reso- 
0^ ' lution simulations" by almost an order of magnitude to examine the convergence of the im- 

^^ . portant physical quantities. The cluster contains ^ 5 million particles within the final virial 

radius, Rmr — 2Mpc (with Hq = 50Kms^^ Mpc^^), and is simulated using a force resolution 
^ . of 1.0 kpc (= 0.05% of Rvir)', the final virial mass is 4.31O"'^^M0, equivalent to a circular ve- 

\J^ ' locity Vcirc = {GAd/Ry^^ ~ 1000 kms^^ at the virial radius. The central density profile has a 

^O ■ logarithmic slope of -1.5, identical to lower resolution studies of the same halo, indicating that 

the profiles measured from simulations of this resolution have converged to the "physical" limit 
down to scales of a few kpc (~ 0.2% of Rvir)- Also the abundance and properties of substructure 
Qv ' are consistent with those derived from lower resolution runs; from small to large galaxy scales 

Q^ . {vcirc > 100 kms"^, m > lO^^AfQ), the circular velocity function and the mass function of sub- 

"^ ' structures can be approximated by power laws with slopes ~ —4 and ^ —2 respectively. At the 

f->i , current resolution, overmerging — a numerical effect that leads to structureless virialized halos 

in low-resolution iV-body simulations — seems to be globally unimportant for substructure halos 
with circular velocities Vdrc > lOOkms"^ (~ 10% of the cluster's Vdrc)- We can identify subha- 
los orbiting in the very central region of the cluster (i?<100kpc) and we can trace most of the 
53 . cluster progenitors from high redshift to the present. The object at the cluster center (the dark 

matter analog of a cD galaxy) is assembled between z — 3 and z — 1 from the merging of a dozen 
halos with Wcirc^300kms~^. Tidal stripping and halo-halo collisions decrease the mean circular 
/\ ' velocity of the substructure halos by « 20% over a 5 billion year period. We use the sample 

of 2000 substructure halos to explore the possibility of biases using galactic tracers in clusters: 
the velocity dispersions of the halos globally agree with the dark matter within <10%, but the 
halos are spatially anti-biased, and in the very central region of the cluster {R/Ryi^. < 0.3), they 
show positive velocity bias {b^ = (Tv3D,haios/o'v3D,DM — 1.2-1.3); however, this effect appears to 
depend on numerical resolution. 

Subject headings: cosmology: theory - dark matter - large-scale structure of the Universe - galaxies: 
clusters - galaxies: halos - methods: numerical 
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^Physics Department, University of Durham, Durham 

City, UK In hierachical cosmological scenarios galaxies 

^Osservatorio di Brera-Merate, Milano, Italy and clusters form in virialized dark matter domi- 

^Astronomy Department, University of Washington, nated halos that are assembled via merging and 

Seattle, USA accretion of smaller structures (White & Rees 

1978, Davis et al 1985; for a recent analysis, e.g. 
Tormen 1997, 1998). Until recently, to what ex- 
tent the subhalos survive within the potential well 



of the larger system has been largely uncertain, be- 
cause cosmological TV-body simulations were not 
able to resolve more than a handful of substruc- 
ture halos (e.g. Carlberg 1994, Summers, Davis & 
Evrard 1995, Frenk et al 1996). Infalling subhalos 
are heated by tidal forces generated by the global 
potential and by mutual encounters and rapidly 
lose a large fraction of their masses; this is a phys- 
ical effect but is greatly enhanced by limited nu- 
merical resolution. The finite resolution sets an 
upper limit to the potential depth of halos - large 
force softening or low numbers of particles per halo 
conspire to produce soft, diffuse substructure ha- 
los that are easily disrupted by tidal forces (Moore, 
Katz & Lake 1996) and lead to structureless viri- 
alised halos. This is the classic overmerging prob- 
lem (White et al 1987). 

Gas physics is not a solution; it is necessary to 
accurately reproduce the dynamics of the (domi- 
nant) dark matter component. It is now clear that 
mass and force resolution of the simulations are 
the key parameters for overcoming the overmerg- 
ing problem (Moore, Katz & Lake 1996, Braincrd, 
Goldberg & Villumscn 1998, Moore et al 1998, 
Ghigna et al 1998, Tormen, Diaferio & Syers 
1998, Klypin, Gottlober, Kravtsov & Kokhlov 
1999a, Okamoto & Habe 1999). Increased reso- 
lution leads to substructure halos {subhalos here- 
after) with higher central densities, enabling them 
to survive. Halos extracted from large cosmolog- 
ical simulations and rc-simulated (see next sec- 
tion) with ^ 10^ particles and force resolution 
<0.01% of the virial radius yield a wealth of sub- 
structure allowing a comparison between the mass 
and light distributions within clusters and galax- 
ies (Ghigna et al 1998, hereafter G98; Okamoto 
& Habe 1999, Moore et al 1999a). Tormen, Di- 
aferio & Syers (1998) have addressed the same is- 
sue using a sample of clusters simulated at lower 
resolution. High-resolution simulations of moder- 
ately large cosmological volumes that retain signif- 
icant amounts of substructure within virialized ha- 
los have also become recently feasible (Klypin et al 
1999a, hereafter KGKK; Klypin, Kravtsov, Valen- 
zuela & Prada 1999b; Colin, Klypin, Kravtsov & 
Khokhlov 1999; Golm, Klypin & Kravtsov 1999; 
see also Kauffman et al 1999a, 1999b and Diafe- 
rio et al 1999). These latter approaches have the 
advantage of providing relatively large samples of 
dark matter halos representing clusters, groups or 



galaxies, but cannot detect systematic biases in- 
troduced by the limited resolution. 

The central density profile of halos is also af- 
fected by numerical resolution. In order to com- 
pare the predictions of hierarchical models with 
observational data on the mass, light and X-ray 
profiles of clusters (e.g. Carlberg et aZ 1996, Carl- 
berg, Yee & EUingson 1997, Small et al 1997, Ab- 
delSalam, Saha & Williams 1998, Adami, Mazure, 
Katgert & Biviano 1998, Allen 1998, Markevitch, 
Vikhlin, Forman & Sarazin 1999), galaxy rotation 
curves (e.g. Moore 1994, Flores & Primack 1994, 
Moore et al 1999b and references therein), giant 
arc properties in gravitationally lensing clusters 
(e.g. Kneib et al 1996; also, Williams, Navarro 
& Bartelmann 1999, Flores, Mailer & Primack 
2000, Meneghetti et al 2000) or constraints on 
processes directly related to the nature of the dark 
matter, such as particle-particle annihilation rates 
(Calcaneo-Roldan & Moore, in preparation), it 
is important to resolve the central structure of 
dark matter halos. Using A^-body simulations 
with 10'' — 10'"' particles per halo for the CDM 
models and its popular variants, Navarro, Frenk 
& White (1996, 1997) found that the profiles of 
isolated relaxed dark matter halos can be well 
described by a "universal" profile (NFW profile) 
from galactic to cluster scales; these results have 
been confirmed by other authors using simula- 
tions of comparable resolution — e.g Cole & Lacey 
1996, Tormen, Bouchet & White 1996. How- 
ever, improving the numerical resolution (Moore 
et al 1998) leads to profiles with central cusps 
significantly steeper than that of an NFW profile 
(p(r) oc r^^ for the latter); halos simulated with 
~ 10® particles have profiles fit by the functional 
form p{r) oc [(r/rs)^-^(l + (r/r^)^'^)]^^ (Moore et 
al 1999b), which has a cusp p{r) ex r^^-^ as r ^ 0. 
(In a recent analysis, Jing & Suto, 2000, find sim- 
ilar results for galaxies and groups, but shallower 
central profiles for a sample of clusters simulated 
in a ACDM cosmology; see our comments in § 4). 

In this paper, we examine how much resolution 
per halo is required to make numerical effects neg- 
ligible for various physical quantities and obtain 
robust results on the halo density profiles and the 
space, mass, velocity distribution of substructures. 
We perform one large (and expensive) simulation 
of a dark matter halo, taking the simulation orig- 
inally analysed by Ghigna et al (1998; hereafter 



G98) and increasing the force and mass resolution 
by almost an order of magnitude. 

The plan of the paper is as follows: In §2, we 
describe the N-body simulations and, in § 3, the 
method used to identify the substructure hales. In 
§ 4, we consider the issue of the typical density pro- 
file of isolated and substructure halos. Section 5 
is devoted to the statisitcal properties of the sub- 
structure, the effects of resolution, evolution and 
environment. We study the distribution of their 
internal velocities and masses, their spatial dis- 
tribution and whether they trace the mass, and 
also the issue of velocity bias. In § 7, we examine 
the relation between the cluster progenitor halos 
at high redshift and the substructure halos at the 
present epoch, and the formation of the central 
"cD" halo. We summarise the results and con- 
clude in § 8. 

2. The A^ body simulations 

Our aim is to achieve very high spatial and mass 
resolution within the dark halo of a cluster drawn 
from a "fair volume" of a standard CDM universe. 
To achieve such resolution, we initially perform 
a simulation of a large volume (50ft.~^Mpc per 
side) of a CDM universe (normalised such that 
(78 = 0.7 and with the shape parameter F = 0.5 
and Hq = 50Kms~^ Mpc~^). The size of this box 
is nearly ten times the turnaround radius allow- 
ing the tidal field and hence infall pattern to be 
modelled correctly. We choose a cluster that is 
virialised by a redshift z — and use a technique 
of "volume renormalization" to obtain higher res- 
olution within the region of interest (see G98 and 
references therein). Our goal is to simulate the 
formation of this cluster such that ^ 10^ parti- 
cles lie within the turnaround radius at the final 
time. (Note that if we simulated the entire volume 
at this resolution we would need to use a total of 
6 X 10^ particles.) Beyond the turnaround region 
the mass resolution is decreased in a series of shells 
allowing us to reproduce the external tidal field 
correctly. The particle distribution is evolved us- 
ing the high performance parallel treecode PKD- 
GRAV, using periodic boundaries and a variable 
timestep scheme based upon the local acceleration 
(Quinn et al 1997, Wadsley, Quinn & Stadel 2000, 
in preparation) . To maintain accurate forces when 
the mass distribution is fairly regular, from a red- 



shift z=69 to z=2 we use an opening angle of 0.4 
and we complete the simulation using an open- 
ing angle of 0.7. We expand the cell moments 
to hexadecapole order. The code uses a spline 
softening length such that the force is completely 
Newtonian at twice our quoted softening lengths. 
The time integration method is an adaptive im- 
plementation of the standard leapfrog integration 
scheme. Individual particle time steps t are cho- 
sen to satisfy the criterium T<ri{e/a)^''^, where e is 
the force softening length of a particle, a the mag- 
nitude of the local acceleration and j] a constant 
determined on the basis of stability and accuracy 
requirements (for these runs rj = 0.2). In order to 
maintain synchronization of the system, the time 
steps are quantized on power-of-two subdivisions 
of the largest time step (the "system" time step 
Tg), i.e. T = Tn = Ts/2" (e.g., Hernquist & Katz 
1989). For these runs, 500-1000 system steps are 
used, equally spaced in time. 

Here we consider two simulations of the cluster 
(the same object as in G98) with largely different 
resolutions, which we label LORES and HIRES. 
In the high resolution region of the box, the par- 
ticle mass is 4.3 x 10^h''^MQ for LORES and 
5.37X lO'^/i-^M© for HIRES; the softening is 1 and 
0.5ft-~^kpc respectively. At the final epoch, the 
cluster has a radius -R200 — l.O/i^^Mpc (i?2oo is 
the radius within which the average density is 200 
times the cosmic density and it is a measure of the 
extent of the virialized region) and a mass within 
such radius M200 = 2.WxlO'^^h-^MQ. With these 
parameters, run HIRES has more than 4 million 
particles within i?2oo , a factor of 8 of improvement 
with respect to LORES. 

The mass growth curve of this cluster is shown 
in Figure 3 of Ghigna et al (1998). The cluster 
does not undergo major mergers since z = 0.5; the 
virial mass increases by <30% since that epoch, 
while the mass within the physical volume of the 
z = cluster (i.e. the mass within R < i?20o|z=o) 
grows by only ^ 10%. The early formation and 
"quiet" recent history of the cluster are conve- 
nient for this study, since they allows us to sin- 
gle out the effects of tides from major accretion or 
merger events. The results presented here are very 
likely typical of observed well relaxed clusters, but 
cannot be easily generalized to objects with more 
turbulent recent histories. 



3. Halo substructure and its identification 

The wealth of substructure that exists within 
the cluster is exemplified by Figure 1. The pan- 
els show maps of the projected density distri- 
bution in a box of side 2i?200 for LORES (up- 
per) and HIRES (lower panel). Each parti- 
cle is plotted using a grey scale according to 
the logarithm of the local density (defined us- 
ing an SPH smoothing kernel over 64 neigh- 
bouring particles using the code SMOOTH of 
Stadel & Quinn 1997, http ref: http://www- 
hpcc.astro.washington.edu/tools). Only regions 
with density contrast (5 > 30 are shown. The 
change of resolution from LORES to HIRES has 
a dramatic effect on the abundance of substruc- 
ture: lower mass objects are resolved throughout 
the cluster and features appear in the very central 
region which was previously smooth. 

The density maps are obtained at redshift z = 
0.5, at which epoch more than 80% of the mass of 
the cluster is in place, yet it is not fully relaxed. In 
Figure 1, we can see three dominant substructure 
clumps (center, north, south-west). Similar fea- 
tures are observed in the light and X-ray distribu- 
tion of observed high-redshift clusters (e.g. Abell 
521; Maurogordato et al 2000). By z = 0, the sec- 
ond largest clump merges with the central object, 
the dark counterpart of a cD galaxy. At z ~ 0.5, 
the cluster's radius and masses are 0.6ft-^^Mpc 
(physical units) and 1.7 x lO^^h^^MQ. 

Identifying the substructure halos is a critical 
step (G98, KGKK). Here we use the group finding 
algorithm SKID, developed by Stadel et al (for 
a full description of the code see http://www- 
hpcc.astro.washington.edu/tools). SKID is a 
"Lagrangian" version of DENMAX (Gelb & 
Bertschinger 1994), the densities being evaluated 
at each particle position using an SPH smoothing 
kernel rather than at the nodes of a grid. The par- 
ticles are moved along the density gradients until 
each oscillates around a point. They are then 
linked using a friends-of-friends algorithm (FOE; 
Davis et al 1985) and the groups (halos) found 
are checked iteratively for self-boundness. 

In G98 we have verified the efficiency of SKID 
and tested the robustness of its estimates of the 
halo structural parameters (see below) against al- 
ternative methods. The ability of SKID to sin- 
gle out particles bound to substructure from the 



diffuse particles bound solely to the entire clus- 
ter (smooth background particles) depends cru- 
cially on the linking length, Ifof, adopted for the 
FOE stage. If IpoF is too small a subhalo will 
be fragmented into smaller units whereas if it is 
too large then halos will be merged together. A 
value Ifof — ^/Srpeak, where Vpeak is the dis- 
tance from the halo's center at which the circular 
velocity profile has its peak Vpeak (see, e.g., G98, 
section 3.3), is usually satisfactory. In our simu- 
lations, the subhalos have a large range of sizes 
and there is not a value for IpoF that is ideal for 
all of them at a time; moreover there are many 
instances of substructure halos that contain their 
own substructure - "halos within halos within ha- 
los". To account for the complete hierarchy of 
substructure present in the current simulation we 
run SKID with three values of Ifof = 1-5, 5 and 
lOZso/t, and combine the outputs avoiding dupli- 
cation. (SKID run with / — lOlgoft provides the 
first level of the "hierarchy" of halos; we then add 
the second level obtained with I = 5lsoft but dis- 
carding those halos that are separated from a first- 
level halo by a distance less than the latter's rpeak ; 
and so forth for the third level. The completeness 
of the final halo catalog in several regions of the 
cluster was verified against the density concentra- 
tions visible in the 3D density map). However, it 
turns out that for the statistics of intermediate- 
to-massive halos (e.g. those with circular veloci- 
ties larger than ^ lOOkms"^) a run of SKID with 
Ifof = ^Qlsoft is adequate (and fast). 

We consider halos with a minimum number of 
member particles of 16 (corresponding to a mass 
of 8.6 • lO^h-^MQ in HIRES). In general, we use 
halos with more than 16 particles when they are 
employed as tracers, but we adopt a minimum 
number of 32 particles when their individual prop- 
erties are relevant. The high resolution region 
analyzed is roughly the final turn-around radius, 
about twice the virial radius. Within this region, 
in HIRES, we identify over 2000 substructure ha- 
los. Within R20Q, there are ^ 1200 halos at the 
final epoch and ^ 1000 at z = 0.5. The inner- 
most halo is at 45ft.^^kpc from the cluster's cen- 
ter although about ten halos lie within a projected 
distance of 45 /i~^kpc. 

We use the output of SKID to determine the 
halo structural parameters: radius rhaio, niass 
Mhaio and "circular velocity" Vdrc (the latter 



is the peak value of the rotation curve v{r) — 
{M{r)/ry^'^, where M{r) is the halo mass within 
a distance r from its center). In § 4.1 of G98, 
we discussed the effects of numerical resolution on 
the completeness of subhalo samples. We mod- 
elled the substructure halos as isothermal spheres 
embedded in a larger isothermal potential — the 
cluster's halo — and used the standard formula 
for tidal stripping by the cluster's potential to es- 
timate the limiting radii (tidal stripping is the 
main cause of subhalo disruption, once the clus- 
ter is in place) . (Using NFW profiles to model the 
halos would not change our estimates significantly; 
see KGKK for a semi-analytical calculation using 
NFW profiles). Applying this simple model to the 
present simulation indicates that the halo samples 
obtained from LORES and HIRES should not be 
affected by major incompleteness for substructure 
halos with Vdrc^ 100 and 60kms^^, respectively, 
assuming pericenters greater than 50/i~^kpc. We 
take these values as quasi-completeness limits of 
the subhalo samples; they mark the values of Vdrc 
below which the samples become apparently in- 
complete. These estimates do not take into ac- 
count halo disruption at redshifts z > 1, before 
the cluster is assembled. For instance, follow- 
ing the cluster's progenitors from high redshift to 
the present shows that the sample derived from 
HIRES may be ^ 20% incomplete starting from 
Vcirc ^ 100 km s^^ (sec § 7). 

4. The density profile of isolated halos 

One of the most fundamental properties of dark 
matter halos is their global density profile. A^- 
body simulations have been extensively used to 
study the structure of cold dark matter halos, 
yielding a variety of important results (e.g. Quinn 
et al 1986, Frenk et al 1988, Dubinsky & Carl- 
berg 1991, Carlberg 1994, Navarro, Frenk & White 
1996, 1997, Cole & Lacey 1996, Tormen, Bouchet 
& White 1996, Fukushige & Makino 1997, Brain- 
erd et al 1998, G98, Tormen, Diaferio & Sycr 1998, 
Huss, Jain & Steinmctz 1999, KGKK, Okamoto & 
Habe 1999). Systematic studies of halos across 
a wide range of masses has revealed a remark- 
able one parameter scaling (NFW profile; Navarro, 
Frenk & White 1996, 1997), with a density profile 
of the form p(r) oc r~^(\ -\- r /rs)~'^ ., where p{r) is 
the spherically-averaged density at the distance r 
from the halo's center and r, a scale radius. The 



halos studied by NFW had w 10^ particles within 
their virial radii and were resolved to a scale of a 
few % of this region. Moore et al (1998, 1999b) 
showed that increasing the resolution by about an 
order of magnitude produces steeper central cusps, 
p(r) oc r~^^^ . Is the current resolution sufficient to 
resolve the central density profile of CDM halos or 
will increasing the resolution make the cusp even 
steeper? 

Figure 2 shows the density profile of the cluster, 
measured in three runs with different resolutions 
lO''-^, lO^"^ and 10^*^ particles within the fi- 
nal virial radius. The lower solid curve shows the 
expected NFW curve for a halo of this mass in this 
cosmology. It provides a good fit to our lowest res- 
olution simulation but underestimates the central 
density of the higher resolution runs. The upper 
solid curve shows the slightly modified profile pro- 
posed by Moore et al (1999; henceforth M99a), 
p{r) oc [{r/rsY'^{l + {f /rsY'^)]~^ , where r^ is a 
scale radius. With increasing resolution, the clus- 
ter's profile appears to converge (asymptotically) 
to M99a's curve. The asymptotic profile is at- 
tained in the regions that have many particles and 
beyond several softening lengths; as it is shown by 
the vertical bars, density profiles measured from 
these A^-body simulations can be "trusted" only 
for scales >6 times the force resolution. This state- 
ment is likely to apply to simulations run with 
different techniques, provided force and mass res- 
olutions arc well balanced, but a detailed study 
of a large number of halos systematically varying 
the numerical parameters, time-stepping criteria 
etc, would be necessary to establish this (see also 
Moore, Katz & Lake 1996 and Moore et al 1998). 

With increasing resolution, the profiles of sub- 
structure halos become closer to those of isolated 
halos. Figure 3 compares the "typical" profile 
of subhalos with Vdrc ^ 200kms~^ in LORES 
and HIRES (obtained for each run by averaging 
the profiles of three subhalos with R ^ 0.5i?2oo 
and Vdrc close to 200 km s^^) with the profile of 
a similar isolated halo, simulated with 10® parti- 
cles (data taken from Moore et al 1999). It is 
interesting to note that the typical subhalo profile 
is steeper in LORES than in HIRES, perhaps as a 
result of numerical relaxation between the hot halo 
particles and cold subhalo particles. This confirms 
that caution is required when studying the prop- 
erties of the profiles of substructure halos even in 



high-resolution simulations (G98; Avila-Reese et 
al 1999, Bullock et al 1999). 

In a recent paper, Jing & Suto (1999) claim 
that the central density profiles of four simulated 
CDM clusters are closer to -1.1 than -1.5 as found 
here. We suspect that this difference may be due 
to the resolution. Jing etal use a Plummer soft- 
ening with quoted resolution 0.005r2oo (their Fig- 
ure 2) and they measure the central slopes using 
the data between 0.007 < r jrmr < 0.02. This re- 
gion may be dominated by dynamics on scales of 
order the softening (as seen above, we find that the 
profiles are reliable at scales that are 3 times the 
scale at which the force becomes Newtonian - our 
spline softening is completely Newtonian at twice 
the quoted values whereas a Plummer softening 
is still "softer" than Newtonian at the resolution 
quoted by Jing & Suto) . It should be noted that a 
steep central profile (with p{r) oc^ ^1-4) was also 
found by Moore et al 1998 using a high resolution 
simulation of a different cluster (the "Coma" clus- 
ter). At the resolutions achieved here and in Jing 
etal, only a handful of systems have been studied 
and the same cluster has not been simulated using 
different codes. We have placed our initial con- 
ditions on our website (www.nbody.net) so that 
other researchers can run this cluster and compare 
with the results presented here and in Ghigna et 
al (1998). 

5. The substructure of dark matter halos 
5.1. Distribution of circular velocities 

The "circular" velocities Vdrc of the subhalos 
are the quantities that can be most easily com- 
pared with observations (although the relation be- 
tween Vcirc and the rotation curve or velocity dis- 
persion of the visible component of a galaxy is 
not straightforward; e.g. Navarro, Frenk & White 
1996). Tidal truncation usually changes the mass 
of a halo by a large amount (often by ^ 80-90%), 
however the peak circular velocity is a fairly stable 
quantity and can be regarded, to some extent, as a 
"label" that can identify a halo for several billion 
years within the cluster environment. 

Major or total halo disruption is significant for 
those few massive halos (with Wcirc^400kms^^) 
that take part in the formation of the cluster 
and undergo similar-mass mergers during the early 
phases of the assemblage, at z ^ 2 (see § 7); 



once the cluster is in place, the dynamical fric- 
tion timescale is much larger than a Hubble time 
even for halos with mass ~ IO^^Mq (for recent 
analyses on dynamical friction see van den Bosch 
et al 1999 and Colpi, Mayer & Governato 1999). 
Most of the substructure halos with small to inter- 
mediate masses survive within the cluster. On av- 
erage, their central masses (the mass within rpeak) 
do not change dramatically and Vdrc varies only as 
the cubed root of the mass. Over 5 billion years, 
Vcirc varies by ~ 20% and mainly for halos that 
spend a large fraction of their orbital period in 
the inner region of the cluster (G98, KGKK; see 
also below). Furthermore, since mergers are rare 
in clusters (G98, Okamoto & Habe 1999), most 
halos conserve their identities. This fact is useful, 
for example, to compare the Vdrc distribution of a 
cluster's substructure halos with that of their pro- 
genitors predicted using the Press-Schcchtcr (PS) 
approximation. 

The differential velocity distribution function 
(VDF), defined as the number of halos per unit 
velocity interval per unit (physical) volume, is 
shown in Figures 4 and 5. The first Figure shows 
the effect of changing resolution. It compares 
the curves obtained for HIRES (solid line) and 
LORES (dotted) at the final epoch (z — 0, up- 
per panel) and for the young cluster (z = 0.5, 
lower panel). We have used all the halos contained 
within the virial radius, i?200, at the two epochs, 
0.975 and 0.6h-^Mpc respectively. The fall off 
at the low-velocity end are caused by incomplete- 
ness due to limited resolution, which dominates 
HIRES and LORES samples for Vdrc < 60kms~^ 
and < lOOkms^^, respectively. LORES's curve 
appears to be slightly affected by resolution up to 
Vdrc ^ 150kms~^; beyond that Vdrc, the results 
for the two runs agree well. In that range, the 
VDF at z = is reasonably well approximated by 
a power-law with exponent —4, motivated by the 
PS approximation for field halos in this mass range 
(e.g. Gross et al 1998, Tormcn 1998, Jenkins et 
al 1999; Klypin et al 1999b). 

In order to single out evolutionary effects, we 
compare directly the halo distributions for run 
HIRES at the two epochs in Figure 5. In the 
main panel of the figure, we compare the VDFs 
measured for the halos within the same physical 
(or proper) volume, as determined by the virial 
radius of the cluster at z = 0. There is virtu- 



ally no difference between the curves for the two 
redshifts (solid and dashed curves) ; the VDF is es- 
sentially "frozen" in time. We find the same result 
for the normalized VDFs, that is, the VDFs ob- 
tained considering only the subhalos within R200 
at each respective epoch and measuring circular 
velocities and distances in units of the cluster's 
circular velocity and virial radius, also at the re- 
spective epoch; the comparison of the two curves 
is shown in the inset of the figure. (For these plots, 
we only show data for Vdrc greater than the quasi- 
completeness limit of HIRES.) 

The amplitude of the VDF of the subhalos 
(those within 0.3-0.5i?2oo) decreases a little over 
time, corresponding to a decrement of Vdrc of 
about 20%. Tidal stripping removes particles or- 
biting at the outskirts of substructure halos but 
tidal heating of the halo cores, from tidal shocks 
near pericenter or from halo-halo encounters may 
also affect their concentrations and accelerate par- 
ticles to the halo boundaries. A 20% decrement 
of Vcirc corresponds to a ~ 50% decrement of 
the mass. This effect may be physical but can 
be enhanced by limited resolution. In fact, the 
LORES simulation exhibits excessive mass loss of 
the inner halos with medium-to-small masses if 
compared with HIRES (see next section). KGKK 
have examined the central mass loss of a model 
subhalo of Vdrc ~ 200kms~^, orbiting in a rich 
cluster with pericenter ~ 150/i~^kpc, varying the 
mass and spatial resolution (see their § 2.4). The 
halo's central mass loss depends on the force soft- 
ening used, if the latter exceeds a certain "op- 
timal" value {^ 1/30 of the halo's tidal radius 
at pericenter, for a Plummcr softening; interest- 
ingly, the dependancc on the particle mass is very 
weak); however, beyond the "optimal" resolution, 
the time evolution of the central mass appears to 
converge to a step function which drops by 20-30% 
at every passage at pericenter, yielding a mass 
loss of ~ 50% over 5 billion years (the orbital 
period in their model is 2 Gyrs). In the HIRES 
simulation, the "optimal" resolution indicated by 
KGKK's test is met e.g. by subhalos of Vdrc ~ 



100 km s 



Sib Tyi 



pert 

200kms-i 



70 h ^kpc and subhalos of 



Vdrc ~ zuUKms " at rperi '^ 50 h ^kpc (we use 



the isothermal model approximation described in 
§ 3 to obtain these estimates) . The fractions of in- 
ner subhalos with Vdrc = 100 and 200kms~^ and 
with estimated pericenters less than 70ft-~^kpc 



and 50/i^^kpc respectively are ^ 20% in both 
cases; therefore the decrease of the VDF from 
z — 0.5 to z = may be affected by numerical 
resolution, but, since the corresponding mass loss 
is close to the value expected from KGKK's test, 
it is probably a physical effect. In any case, it sets 
an upper limit. (For comparison, in LORES, the 
"optimal" resolution for rp^ri ^ 70/i^^kpc is met 
only if Wcirc^200kms^^). 

The shape of the VDF does not depend signifi- 
cantly on the environment. This is also shown by 
Figure 5, where we plot the VDF for the inner sub- 
halos {R/R20V1 < 0.5; upper dotted line) and those 
in the cluster's periphery (1 < R/R200 < 2.5; 
lower dotted line) . A power-law dN/dvdrc oc v^^^^ 
is always a reasonable fit for fci,c^l00kms^^. 

5.2. Mass distribution function 

Knowledge of the mass distribution function 
(MDF) is useful to estimate the disruptive power 
of halo-halo collisions in dense environments. For 
example, modelling this process is necessary to un- 
derstand the origin of intracluster light and the 
morphological evolution of galaxies in clusters and 
groups ("galaxy harassment"; Moore et al 1996). 
The mass and radii of substructure halos can 
also be probed directly using mass-reconstruction 
methods based on observations of arcs and arclets 
in clusters (Natarajan et al 1998). 

The MDF provides information independent of 
that derived from the VDF. The mass of a halo of 
given Vdrc depends on its orbital history, which de- 
termines the extent of tidal disruption that it has 
suffered. This strongly correlates with the smallest 
pericentric distance reached by a halo (G98) and, 
to a lesser degree, with the number of orbits com- 
pleted (since a halo loses mass very quickly the 
first time it passes through pericenter and more 
slowly every time it completes an orbit; KGKK). 
However, since halos of different masses but with 
the same pericenter lose similar fractions of mass 
and the distribution of pericenters does not de- 
pend significantly on mass (G98), we still have 
roughly Mhaio oc Wcirc; even for the subhalos. 

Figure 6 plots the MDF for the substructure 
halos (defined as the number of halos per unit log- 
arithmic mass interval per unit physical volume). 
Results are shown for both runs at the two red- 
shifts considered previously. At z = 0.5, the agree- 



ment between HIRES and LORES is good. The 
mass limit for which both runs are close to com- 
pleteness (in the circular velocity distribution) is 
>5.6-lO^°/i~^M0, the typical value for a halo with 
Vc„c = lOOkms-1 (see Figure 20 of G98). At 
X = 0, the MDF for the LORES cluster is ~ 50% 
lower than for the HIRES simulations. This dif- 
ference, as for the VDF, is caused by numerically 
enhanced mass loss near the cluster's center; it 
disappears if we exclude the halos in the central 
region, i?/i?200 < 1/2, from the counts (Figure 7). 
As discussed in the previous section, our HIRES 
run is not affected significantly by artificial mass 
loss for Vcirc>100kms-\ i.e M>6 • IG'^^H-'^Mq. 

As found for the VDF, the MDF does not ex- 
hibit significant evolution within a physical vol- 
ume during the lifetime of the cluster (Figure 8, 
main panel). At both redshifts, the MDF is 
close to a power-law dn{M)/dM ex Af^^, for 
M/ja;o^lO^^/i~^M0, but becomes shallower for 
lower masses, dn{M)/dM oc M~^-'' (in the fig- 
ure, we plot dn/dLogioM). The main panel of 
Figure 8 also shows the MDF for inner and outer 
halos (symbols as for Figure 5); all the curves 
are similar, with slopes varying between —0.8 and 
—0.9. The inset shows the normalized MDFs for 
the cluster subhalos at the two redshifts, obtained 
in a similar way as the normalized VDFs of Fig- 
ure 5 in the previous section. 

The mass bound to substructure halos is a lit- 
tle more than 10% of the total cluster mass, with 
more fractional mass in halos at early rather than 
at late epochs. The value varies very slowly with 
resolution, but it can oscillate significantly since 
a large fraction of the mass belongs to the two 
or three most massive subhalos (vdrc ^ 400 — 
500 kms^^) that contribute « 5%. If they are 
close to the cluster's center or if they possess their 
own substructure, the masses measured by SKID 
are quite sensitive to the linking length used and 
can differ significantly from estimates obtained 
adopting the spherical overdensity method (see 
G98), sometimes by a factor of 2. For instance, us- 
ing the spherical overdensity method for LORES 
we find a global fraction of mass in subhalos of 
~ 13% (G98), while using SKID's masses this 
value reduces to ~ 11%. If we exclude the largest 
substructure halos, different mass estimates yield 
halo mass fractions differing only by about 1%. 

Figure 9 plots the fraction of cluster mass that 



is in subhalos with circular velocity above a given 
value Vcirc, but below 400kms~^ The difference 
between the two runs never exceeds 1%. Also, 
between z ~ 0.5 and z — 0, the mass fraction 
decreases by only ~ 1% (note that the mass frac- 
tion is computed relative to the virial mass, M200J 
at the respective epoch, namely 1.7 and 2.15 x 
10^'^h~^MQ). Where the halo samples are close 
to completeness, the results are well represented 
by the function Mhaios{vcirc > Vcirc)/M2ao - 
-0.091ogio[wcirc/(kms-i)] + /?, with /3 = 0.23- 
0.22. The mass fraction varies roughly linearly 
with clusterccntric distance (G98, Figure 11). 

5.3. The spatial distribution of substruc- 
ture halos 

In this section, we examine the bias existing be- 
tween the spatial distribution of particles (i.e., the 
mass) and that of the subhalos. Here, the bias 
is defined as the ratio b{R) = nhai(,R) / npart{R) , 
where nhai{R) and UpartiR) are the number den- 
sity profiles of subhalos and particles as functions 
of clusterccntric distance R. For convenience we 
normalize all quantities to the cluster parameters, 
i.e. we measure lengths in units of the virial ra- 
dius, -R200J and Nfiai and Npart as fractions of 
the total numbers of subhalos and particles within 
i?200- The halo distribution will be (positively) bi- 
ased \ib{R) increases as R decreases, anti-biased in 
the opposite case. In G98 we found that the halos 
are anti-biased; Colm et al (1999) also find anti- 
bias for dark matter halos in large volume simu- 
lations of four cosmological models (they measure 
the bias defined by the ratio between the 2-point 
correlation functions of mass and dark matter ha- 
los, as it is customary for large scale structure 
studies). To study this problem it is particularly 
important to pay attention to the limitations im- 
posed by the finite resolution, since overmerging 
— preferentially erasing subhalos in the central re- 
gions of the cluster — always introduces anti-bias 
in the subhalo distribution (or enhances it). 

We compare the normalized (or fractional) 
number density profiles of halos and particles as a 
function of clusterccntric distance R in Figure 10, 
for the redshifts z ~ 0.5 and z = (top and 
bottom panel respectively). We plot the num- 



bers of halos with ' 



>v. 



LIM 



in spherical shells 



of radius R divided by the volume of the shell 

(in units of -R20o)j divided by the total number 



(-/V200) of cluster subhalos above the circular ve- 
locity limit (the "cluster" subhalos are those con- 
tained in i?200 at the respective epoch; this yields 
188 and 219 subhalos if w^^^f = 80kms-\ and 
104 and 110 if vj:l^ = 100 kms-^ for HIRES, 
while for LORES the corresponding numbers are 
165 and 158, and 116 and 100). For the mass, we 
plot the particle number density profile divided 
by the number of particles within i?200- In each 
panel, the long dashed curves are the normalized 
mass profiles; the two other curves are for the ha- 
los of HIRES {solid with errobars) and LORES 
(dotted). For z ^ 0.5, we set v^lJ^ = 80kms~\ 
which yield good statistics; for z = we use 
v^lJf = 100 km s^^, which has similarly good 
statistics. Varying the value of v^l^ does not 
make any difference, as long as there are enough 
halos to obtain significant measures. 

Within the cluster, at z — 0.5, the halos have a 
number density profile similar to that of the mass, 
but the latter profile declines more steeply than 
the former approaching the virial radius and be- 
yond (i.e., b{R < i?20o) < b{R > -^200)- This 
means that the population of halos (for a circu- 
lar velocity limited sample) in the volume encom- 
passing the cluster is globally anti-biased - with 
proportionally fewer halos where the particles are 
more clustered. At z = 0, there is also anti-bias, 
the mass being clearly more concentrated than the 
halos; for R ^ 10% of i?2oo, there is a factor of 2 
of difference between the two profiles (and is even 
larger if we require that the two profiles overlap 
beyond i?20o)- At both redshifts, there are no sig- 
nificant differences between the curves for HIRES 
and LORES. Therefore the halo number density 
profiles measured from this simulation should not 
be significantly affected by residual overmerging 
(see also § 6). 

Figure 11 compares directly the halo number 
density profiles a,t z ~ 0.5 and z = in the phys- 
ical (proper) volume centered on the cluster (for 
both samples, we set v^/^^ = 80kms~^). The ha- 
los are clearly more concentrated in the volume 
of the newly formed cluster. Infalling material 
subsequently enhances the halo number density in 
the outer region, i?>600/i~^kpc (the virial radius 
of the cluster at z = 0.5), while tidal disruption 
in the inner region reduces the number density of 
halos with circular velocities above the limit. It 
is interesting to note that allowing for a ^ 20% 



decrease of a halo's Vdrc between z = 0.5 and 
z — (because of central mass loss, as observed in 
the VDF analysis in § 5.3) is sufficient to recover 
most of the central steepness of the z = 0.5 pro- 
file; this is shown by the dotted line in the Figure. 
This indicates that, once the cluster is in place, 
tidal mass loss alone (instead of total halo evapo- 
ration) causes the measured decrement of the sub- 
halo number density towards the cluster's center. 
Apart from tidal mass loss, mergers between 
halos during the cluster assembly (at z > 0.5) are 
probably a major cause of the anti-bias between 
halos and particles. At z = 0.5, the uniform anti- 
bias of the cluster subhalos (with respect to the 
'field'; as seen in Figure 10, upper panel) may re- 
flect the fact that the cluster has only recently 
being assembled through the merging of group- 
sized halos of similar masses, so that mergers and 
tides have affected its subhalo population more 
uniformly. Note that the galaxies that would re- 
side in the subhalos could have a different bias, if 
their circular velocities are not directly affected by 
the central mass loss of their dark halos. As a final 
word of caution, we point out that the issue of bias 
would be better addressed with a statistical sam- 
ple of clusters and a more representative sample 
of field halos. This approach requires a trade-off 
between statistics and resolution but considerable 
progress has recently been made (Colin et al 1999; 
see also Kauffmann et al 1999a,b, Diaferio et al 
1999, Benson et al 2000). 

5.4. Velocities of halos and dark particles 

The issue of velocity bias — the difference be- 
tween galaxy and dark matter (DM) velocities — 
was first addressed by Carlberg & Couchmann 
(1989) as a means of reconciling low measures of 
fi from the observed velocities of galaxies in clus- 
ters. There are two types of velocity bias consid- 
ered in the literature, the one-point velocity bias 
{by = (Jv.haios/'^v.DM) which compares the veloc- 
ity dispersions of galaxies and dark matter par- 
ticles (e.g. Carlberg & Dubinski 1991, Carlberg 
1994) and the two-point or pairwise velocity bias 
(6„^i2) which compares the relative velocity disper- 
sions in pairs of objects (e.g. Couchmann & Carl- 
berg 1992). Here we consider the first type. A re- 
cent paper (Colin, Klypin & Kravtsov 1999; CKK 
hereafter) focussed on the different forms of veloc- 
ity bias for numerical simulations of the ACDM 



cosmology. The situation is still rather confusing. 
CKK find that the dark matter halos within clus- 
ters are positively biased, with by ~ 1.2-1.3 (simi- 
lar results were also found by Diaferio et al 1998 
and Okamoto & Habe 1999), while KGKK find 
bv^^ for a similar set of data (see their Figure 15). 
In G98, we did not detect velocity bias. The dif- 
ference may be partly due to the fact that G98 
and KGKK measure the global b^ using the whole 
sample of cluster subhalos, while CKK consider 
the radial profile bv{R). We re-examine the issue 
here in more detail. Note that the comparison 
between the spatial distribution of the halos and 
that of the mass considered in the previous sec- 
tion, does not provide sufficient information, since 
the scaled profiles of differently biased tracers may 
not differ much (Carlberg 1994). 

We compare the 3D velocity dispersion {(Tv.sd) 
profiles of halos and dark matter for run HIRES 
in Figure 12 (the averages are taken in spherical 
shells of width 250ft.~^kpc; at the center, (Tv.sd 
for the DM drops sharply from ^ 1300 km s^^ at 
R ~ 50ft.-ikpc to lOOOkms-i at i? ~ lO/i-^kpc, 
if radial bins of width 10 /i^^kpc are used, as it is 
shown by the light solid curve in the lower panel of 
the figure). We consider the cluster at z = (lower 
panel) and also the output at z = 0.1, separated 
in time from the final epoch by less than one typ- 
ical subhalo orbital period; we also subdivide the 
halos in three subsamples according to the values 
of their circular velocities (as listed in the figure) . 
In this way, we should be able to single out signifi- 
cant differences between the (3D) velocities of ha- 
los and dark matter from statistical flukes, partly 
overcoming the drawback of dealing with only one 
cluster. The samples with Vdrc > 150kms~^ have 
the worst statistics, especially for R > -R200 where 
they provide ^ 5 halos per bin. 

With quite large oscillations, the (Jv.sd of the 
halos is consistent with that of the dark matter 
for i?>300/i^^kpc. Within the cluster, the data 
show a significant (positive) bias in the innermost 
bin, with by — 1.2-1.3. The maximum amplitude 
of bias in our data is similar to that reported by 
CKK but the signal is globally weaker (CKK have 
by ~ 1.2 at R/R2oa ^ 0.7, that is >600/i"^kpc for 
the cluster studied here). We also show the veloc- 
ity dispersion profiles of halos and particles for run 
LORES, in Figure 13. The trend of increasing by 
with decreasing R is more apparent in this case. 



Here, the sample with the lowest circular veloci- 
ties is largely incomplete due to limited resolution; 
interestingly, it exhibits a more marked bias than 
the other samples. It is also interesting to note 
that there are no obvious biases in the distribu- 
tion of the orbital parameters of the halos and the 
particles (we have computed the orbits approxi- 
mating the cluster potential to a static spherical 
potential obtained from its density profile, as in 
G98, § 4.6). 

When we average over the virial volume of the 
cluster, the velocity bias of the halos in the cen- 
tral region has a small impact; in the three cases 
considered previously, the global bias within i?200 
is always less than 1.1 (although it remains posi- 
tive) . Also considering only massive halos {vdrc > 
200 kms^^) does not yield any significant bias. In 
this case, anti-bias would be expected if dynamical 
friction were efficient in slowing down such halos, 
but the friction time for subhalos within a clus- 
ter is almost always larger than the Hubble time 
(see Colpi, Mayer & Governato 1999 for a recent 
analysis). Even restricting the samples to half the 
virial radius does not make a significant difference 
(in one case, for HIRES at z ~ 0.1, there is anti- 
bias, but the other two yield positive or no bias, 
unless the "cD=global potential" is included). 

In conclusion, we find evidence of positive ve- 
locity bias in the central region of the cluster but 
the signal is much weaker than CKK's detection. 
It seems that limited resolution enhances the sig- 
nal; overmerging accentuates the impact of phys- 
ical processes like mergers and total halo disrup- 
tion that could lead to positive velocity bias. A 
larger set of high-resolution simulations would be 
welcome to address these issues with better statis- 
tics. 

6. Survival of substructure halos 

As seen in § 5, over the lifetime of the cluster 
the abundance of substructure does not appear 
to change significantly; however substructure is 
erased also at earlier redshifts. To examine this 
issue, we have identified the dark matter halos that 
are progenitors of the cluster at z = 1 and z = 3 
(using the public friends-of-friends code available 
from http://www-hpcc.astro.washington.edu/tools, 
run with a linking length of 0.2 times the mean 
interparticle separation). We define as progeni- 
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Fig. 1. — Maps of the cluster's density as seen in 
run LORES (upper panel) and HIRES, which has 
8 times better mass resolution. The change in the 
appearance between the two runs is mainly due 
to HIRES 's ability of resolving further down the 
substructure mass function. 



Vcirc, min 


Nhalos 


Ptraced 


PcD 


230kms-i 
lOOkms-1 
50kms-i 


15 
120 
550 


90% 
100% 
82% 


-10% 

0% 
0% 



Table 1: Results of tracing the cluster's progenitor 
halos at 2 = 1 to the substructure halos at the final 
time {descendents) , for three ranges of mass of the 
progenitors (corresponding to more than 10, 000, 
between 10, 000 and 1, 000, and between 1, 000 and 
100 particles). The first and second columns re- 
port the minimum circular velocity of the progen- 
itors in each mass range and their number; the 
third and fourth columns give, respectively, the 
percentage of progenitors in each range with a de- 
scendent at z = other than the "cD" object, and 
the percentage of progenitors that merge with the 
"cD" . (These values are obtained requiring for the 
descendents a minimum number of traced particles 



N^ 



equal to 1% of the particles in the progen- 



itor and not less than 4; there is no difference if 
we set Np^rnin to be at least 16). 



Fig. 2. — The density profile of the cluster mea- 
sured in three runs with increasing resolutions 
(from triangles to squares to circles). In the best 
run, the cluster contains over 4 million particles 
and the force resolution is 0.05% of the cluster's 
virial radius. The curves are an NFW profile 
(lower curve) and a fit with the profile of Moore et 



-i.b\ 



at 



al (1999a), which rises more steeply (oc r 
the center than the NFW profile (oc r'^). With 
increasing resolution, the cluster's profile contin- 
ues to approach M99a's curve i.e. this appears to 
be the asymptotic profile in the limit of infinite 
resolution. The vertical bars mark the radii at 
which the measured profiles are no longer affected 
by finite numerical resolution. 



'^circ, min 


-^^ halos 


Ptraced 


PcD 


300kms-i 
150kms-i 
70kms-i 


23 

180 
1070 


40% (40) 
90% (90) 
70% (61) 


60% (60) 
10% (7) 
3% (0.1) 



Table 2: Results of tracing the cluster's progenitor 
halos at z == 3 to the substructure halos at the 
final time. The mass ranges and the quantities in 
the columns are the same as in table 1 (in the last 
two columns, the values in brackets give the per- 
centages obtained requiring Np,min to be at least 
16). 



11 
















: 


10» 


_ 


^^ 






_ 


105 


r 








-= 


10" 








\ 


1 


000 


^ - 


Subhalo (Lores) 

Subhalo (Hires) 


\ 


-. 


100 




Is 


olated halo (v^ = 


220kms-') \ 


^ 







0.01 



0,1 



R/R„ 



Fig. 3. — The typical density profile of a substruc- 
ture halo with Vdrc ~ 200knis~^ (and cluster- 
centric distance '^ 0.5i?2oo) in the two runs, com- 
pared with the profile of an isolated halo (simu- 
lated with much higher resolution) . Increasing the 
resolution brings the profiles of substructure halos 
closer to those of their isolated counterparts. 
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Fig. 4. — Velocity distribution function (VDF) 
of the substructure halos, at two redshifts (lower 
panel: the final epoch, z = 0; upper panel: the 
young cluster at z = 0.5, when it is ~ 1 billion 
years old). The figure shows the effect of increas- 
ing the resolution by a factor of 8 in mass (from 
LORES, dashed curve, to HIRES, solid curve). 
The errorbars represent (1-cr) Poisson errors on 
the counts. The curves agree well where we expect 
both runs to be close to completeness (roughly 
Vcirc > lOOkms^^). The fall off at the low ve- 
locity end is caused by finite numerical resolution. 



Fig. 5. — The VDF for the cluster's substruc- 
ture halos at z = (solid line) compared with 
that of the halos a,t z — 0.5 contained within 
the same physical volume (dashed), using in both 
cases i?20o|z=o £^s limiting distance (main panel). 
The inset shows the VDFs obtained considering 
only the halos within the cluster's virial radius at 
z — (dashed) and z = 0.5 (dotted) and and mea- 
suring Vcirc in units of the cluster's "circular veloc- 
ity" at the respective epoch. In both cases, there 
are no significant changes in the shape and ampli- 
tude of the VDF during the lifetime of the cluster. 
The two dotted curves in the main panel show the 
effect of changing environment inside and around 
the evolved cluster (z = 0); they are for inner sub- 
halos (i?/i?200 < 0.5; upper curve) and "periph- 
eral" halos (1 < R/R200 < 2.5; lower). Within 
and around the cluster the shape of the VDF is 
very similar. 
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Mass Distribution Function (MDF) 




Fig. 6. — The mass distribution function (MDF) 
of substructure halos at two redhifts (lower panel: 
2; = 0; upper panel: z — 0.5) for the two runs 
(LORES, dashed curve; HIRES, solid curve). The 
errorbars are (l-c) Poisson errors on the counts. 
As for the VDF, at z — 0.5, the curves agree well 
where we expect both runs to be close to com- 
pleteness (M>5.6 • W^^H-'^Mq). At z = 0, the 
deficiency of LORES 's MDF is due to numerically 
enhanced mass loss of medium to small mass halos 
in the central region (see Figure 7). 
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MDF, dependence on time and environment 




Fig. 8.— Main panel: the MDF a,t z = for the 
cluster's halos (solid line, HIRES data) compared 
with that of the halos at z = 0.5 contained in 
the same physical volume (dashed line; the curves 
are shown only in the range where they are not 
significantly affected by resolution). The mass 
function is "frozen in time" and its shape does 
not depend much on the environment either, as it 
can be seen by considering only the inner subhalos 
(/J/-R200 < 0.5; upper dotted curve, main panel) 
or those in the "periphery" (1 < i?/i?200 < 2.5; 
lower dotted curve, main panel). The inset shows 
the normalized MDFs at the two redshifts, ob- 
tained in a similar way as the normalized VDFs 
of Figure 5. 
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Fig. 7. — Same as in the lower panel of Figure 6, 
but here the halo samples have been constructed 
requiring R > l/2i?2oo- The curves for HIRES 
(solid) and LORES (dashed) are close in this case. 
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Fractional number density vs R (z = 0.5) 
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Fig. 9. — Fraction of cluster mass bound to subha- 
los with circular velocities exceeding a given value 
Vcirc, at z = 0.5 (upper panel) and z = (lower), 
for HIRES (solid line) and LORES (dashed). We 
also set an upper limit of 400kms~^ to exclude 
the most massive halos whose masses have a rel- 
atively large scatter; they contribute alone an ad- 
ditional ^ 5%. The light dotted line is a fit to the 
z = 0.5 curves for Vdrc > lOOkms"^ (the flatten- 
ing of the curves at the low- velocity end is due to 
the incompleteness of the samples). The bound 
mass fraction is oc log (ucirc), thus varies slowly 
with resolution. 
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Fig. 10. — Fractional number density of halos and 
particles as a function of clustercentric distance 
R, a,t z = 0.5 (upper panel) and z = (lower 
panel); lengths are in units of the virial radius 
i?2oo at the respective epoch; A^20o is the number 
of halos or particles within i?200- In the top panel, 
the solid curve (with 1-a Poisson errobars) is for 
HIRES's subhalos with Vdrc > 80kms~^; the dot- 
ted curve is LORES's equivalent. The dashed line 
is the particle (i.e. the mass) profile. The subha- 
los are always antibiased with respect to the mass, 
although at z = 0.5 their profile is similar to that 
of the mass within the cluster. 
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tors those halos that contribute particles to the 
virial volume of the cluster at the final epoch. We 
have then "traced" them to the substructure halos 
identified at z = 0, by comparing the indexes of 
their particles. For the z ~ subhalos we use only 
the 'core' particles, i.e. those that are separated 
from the halo centers by a distance less than the 
value of Tpeak " this (conservative) approach avoids 
ambiguities when there is a hierarchy of substruc- 
ture. For the central object, we adopt a radius 
rcD — 25/i~^kpc, a value that should single out 
the progenitor halos that merge to form the "cD" 
core against those that only supply particles to its 
halo through tidal stripping. 

We accept descendents only if they have a num- 
ber of traced particles larger than 1% of the par- 
ticles of the progenitors and not less than a min- 
imum number, Np^im which we varied between 
4 and 16 (the number of particles traced can be 
small for small halos since we only consider their 
'core' particles). In our approach, a progenitor 
with no descendent subhalos at z = should have 
been completely disrupted with its particles lost 
in the diffuse particle background of the cluster 
or, possibly, in the outer parts of large subhalos 
(visual inspection of 10 cases confirms this). Such 
halos may survive if the numerical resolution were 
increased further, so that the fraction of progen- 
itors with no descendents is an estimate of the 
importance of residual numerical effects. Cases in 
which two or more progenitors have a common de- 
scendent, i.e. the progenitors' particles contribute 
to the desccndent's core, point to true mergers. 

Table 1 reports the results of the tracing for the 
z = 1 progenitors in three mass ranges. Table 2 
is the analog for z = 3. In the first column, we 
show the minimum value of the circular velocity 
for the progenitors in each mass range. All the 
z = 1 progenitors with Vdrc > 100 kms^^ have a 
descendent at z = 0; the tracing is still essentially 
complete at Wcirc^70kms~^ {ptraced = 98%) and 
deteriorates only for Vdrc approaching 50kms~^ 
{ptraced — 80%). The Central object accretes only 
another halo {oivdrc ^ 300kms~^). These results 
are stable against varying Nmin ■ The loss of halos 
from z ~ 1 to the present appears negligible for 

Vcirc > 100 km S^^. 

The most prominent feature of the tracing of 
the z = 3 progenitors is the large fraction of 
massive halos that merge with the central object: 



about 25 progenitors with Vdrc > 150, a dozen 
for Vdrc > 300kms~^. There is virtually no con- 
tribution to the central object's (core) mass from 
progenitors with Vdrc < 150. Essentially all the 
progenitors with Vdrc > 150kms~^ have a part- 
ner at z = 0; the tracing starts deteriorating at 
Wcirc^lOOkms"^ where the percentage of "child- 
less" progenitors is 18%, independent of Nmin in 
the range 4-16. Successful tracing with the "cD" 
might still hide overmerging, but it seems safe to 
consider the mergers of the massive halos as real, 
since they occur before the cluster formation and 
involve objects of comparable masses (for which 
the dynamical friction time is small). Also the 
distribution of their particles in the z = cluster 
is well concentrated within the cD. We have ex- 
amined the distribution at the final epoch of the 
particles in the cores of the z = 3 progenitors with 
Vdrc > 200kms~^ that merge with the cD (we use 
rcore — 20kpc comoving); about 80% of such par- 
ticles end up within 50ft-~^kpc from the cluster's 
center at z = 0. Very likely further increases of 
resolution, and the inclusion of baryons with cool- 
ing, will reveal substructure within a cD's radius; 
e.g. the substructure already present within the 
progenitors at z = 3 could survive and be found 
within the cD. However, these subhalos would lose 
a very large fraction of their masses and should 
have small circular velocities and low mass-to-light 
ratios. 

In conclusion, overmerging effects in this sim- 
ulation appear negligible for large and intermedi- 
ate halos; loss of substructure may still affect ha- 
los starting from Vdrc ^ 100 kms^^ at the ^20% 
level. This does not seem to be important for 
most of the statistical properties of the substruc- 
ture considered in the previous section, in view of 
the strict similarities of the results for LORES and 
HIRES (except possibly for the velocity bias). 

7. Summary and Conclusions 

We use "N-body" simulations to follow the col- 
lisionless evolution of the dark matter component 
of a rich galaxy cluster. Increasing the force and 
mass resolution (Figure 1) by an order of magni- 
tude has a dramatic effect on the abundance of 
substructure and allows us to assess the biases in- 
troduced by low resolution studies. We find the 
following main results: 
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• We compare the density profile of the clus- 
ter halo simulated with 2 x lO**, 6 x 10^ 
and 4 x 10^ particles within the virial ra- 
dius. At our highest resolution the force 
softening is 0.5/i~^kpc (0.05% of the virial 
radius) and the profile is well fit by p{r) oc 
[(r/r^)i-''^(l + (r/rs)i-5)]-i. This agrees with 
the results of Moore et al (1999b) for lower 
resolution simulations, and suggests that we 
may have converged on the asymptotic cen- 
tral slope. The best fit NFW profile to the 
data has residuals of order 30%. (In a recent 
paper, Jing & Suto, 2000, find shallower cen- 
tral cusps {p(r) oc r^^'^) for four high reso- 
lution simulations of clusters in the ACDM 
model; the reasons of the disagreement are 
not yet clear, see § 4). 

• For the first time we can confidently re- 
solve the cluster density profile beyond 
« Q.2%Ryir = 2/i"ikpc. Between 2/i"ikpc 
and 40 /i~^kpc the density profile has a cen- 
tral cusp with a slope p{r) oc r^^-^ ± 0.05. 
This disagrees with the central mass profile 
of CL0024+1654 recovered from a strong 
Icnsing analysis which indicate a constant 
density core over the same radial range 
(Tyson et al 1998). 

• The density profiles of substructure halos 
has not yet "converged" . Lower resolution 
studies found that halos within clusters had 
steeper density profiles than field halos - this 
may be due to particle-halo heating of sub- 
structure halos. 

• The distribution function of substructure 
circular velocities (VDF) is essentially in- 
variant over the lifetime of the cluster (^^ 5 
billion years) and is quite close to a power- 
law dn{vcirc) / dvcirc OC u^'*, with little de- 
pendence on environment (§ 5.1). Clusters 
simulated with ^ 10^ particles can reliably 
measure this function for halos with circular 
velocity above lOOkms"^. (Note that the 
cluster forms early, at z>0.5, with no major 
mergers since; the evolution of the properties 
of the substructure may be different in ob- 
jects with more turbulent recent histories.) 

• The distribution of halo masses (MDF) 
shows the same "invariance properties" as 



the VDF -virtually no evolution and little 
dependence on environment- and is close to 
a power-law dn/dM oc M~^ for halo masses 
>lO"/i-iM0 (a fraction 5 x lO"'^ of the 
cluster's virial mass); tidally driven mass 
loss changes the MDF self-similarly (§ 5.2). 
Clusters simulated with ~ 10^ particles can 
reliably measure this function for halos with 
masses above 6 x 10^^h~^MQ. 

• The substructure halos of circular veloc- 
ity limited samples are spatially anti-biased 
with respect to the underlying mass distri- 
bution (see § 5.3 for definitions). At an early 
epoch (z = 0.5), the anti-bias is uniform 
within the cluster (halo and particles have 
number density profiles of similar shape); 
at the final time, the anti-bias of the ha- 
los increases approaching the cluster's cen- 
ter. Halo-halo collisions and tidal mass loss 
(which decrease a halo's Vdrc, on average by 
~ 20% over 5 billion years) are the likely 
source of the anti-bias (§ 5.3). 

• The 3D velocity dispersions of halos and 
dark matter particles are generally con- 
sistent within and around the cluster, ex- 
cept in the central region {R < 0.3i?2oo — 
300/i^^kpc) where the halos are "hotter" 
than the particles with a (positive) veloc- 
ity bias by = ay^haios/(^v,DM = 1.2-1.3). 
The sign and maximum amplitude of the 
bias is similar to that reported by Colin et 
al (1999), but the radial range over which 
6t, > 1 is much more limited (their detection 
extends to i? ~ 0.8i?2oo)- The magnitude of 
the bias appears to be enhanced by limited 
resolution since "overmerging" erases prefer- 
entially the central cluster halos. When we 
average over -R200: the velocity dispersions 
of halos and particles do not differ by more 
than 10%. 

• Following the cluster progenitors from high 
redshift to the present shows that, at the 
current resolution, overmerging should be 
globally unimportant for subhalos with 
Wcirc^lOOkms^^, although a fraction of 
high-redshift halos may be lost for Vdrc 
close to 100 km s^^. All the cluster pro- 
genitors with Vcirc<^70]<iiRS~^ at z = 1 have 
a descendant among the subhalos at z = 0; 
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about 20% of the progenitors a,t z — 3 with 
Vcirc ~ lOOkms"^ may be artificiaUy dis- 
rupted. The "cD" object at the cluster's 
center is assembled at high redshift (z ^ 3 
to 1) through the merging/ accretion of a 
dozen halos with Wcirc^300kms~^; it also 
accretes one large halo within the cluster 
(at z < 0.5). 

It is now clear that the hierarchical model 
for structure formation can naturally form mas- 
sive dark matter halos that contain a wealth of 
substructure resembling observed galaxy clusters; 
high resolution numerical simulations allow us to 
obtain robust statistical measures of the proper- 
ties of the substructure and make detailed predic- 
tions. Since substructure halos are expected to 
host galaxies in a cluster and galaxy satellites in 
a galaxy, the comparison with the observed dis- 
tribution of those objects can produce important 
tests for the cosmological models. For example, 
whereas the MDF obtained from cluster simula- 
tions in the CDM universe easily reproduces that 
derived from the observed circular velocity func- 
tion of cluster galaxies (Moore et al 1999a), the 
predicted abundance of substructure of galactic 
halos is not matched by the observed population 
of satellites of the Milky Way and Andromeda 
(Klypin et al 1999b, Moore et al 1999a). 
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Halo number den 




Fig. 11. — The number density profiles in physical 
(proper) volume using halos with Vdrc > "Vcirc = 
80kms^^ at z = (solid line) and z — 0.5 
(dashed) . The halo number density decreases in 
the central region (-R<i?20o/3). The dotted line 
shows the effect of having a variable v^^^^ that 
drops by ^ 20% for R < -R200/4 (precisely, by 
15% in the radial range 0.125 < R/R200 < 0.25 
and by 25% for i?/i?200 < 0.125). 



This 2-coluinn preprint was prepared with the AAS lATJiiX 
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Fig. 13. — Same as in Figure 12 but using data 
from LORES at z = 0. The positive bias of the 
halos in the central region is more marked than 
for HIRES. 



Fig. 12. — 3D velocity dispersion profiles for the 
dark matter (solid curve) and for three samples 
of halos (the points on the dotted curves) selected 
according to their circular velocities as outlined 
on the figure {vdrc in kms^^). The data are for 
HIRES at the redshifts z = 0.1 and z = 0; R200 
is respectively 0.88/i~iMpc and 0.975 /j-^Mpc. 
(The curve for Vdrc > 150 does not extend to the 
innermost bin, since the latter contains only two 
such halos - the "cD" has been excluded (it is al- 
ways at rest since it simply the central smooth 
cluster potential; the profile for the dark matter 
using fine bins near the cluster's center is shown by 
the light solid line in the lower panel). The halos 
in the central region (i?<300 h^^kpc) are "hotter" 
than the DM, with 6^,-1.2-1.3. 
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